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We present a study of the achievable rates for classical information transmission via a lossy bosonic channel 
with memory, using homodyne detection. A comparison with the memoryless case shows that the presence of 
memory enhances the bit rate if information is encoded in collective states, i.e. states which are entangled over 
different uses of the channel. 
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Q\ . I. INTRODUCTION 

O ; 

■ One of the tasks of quantum information theory is to evaluate the capacity of quantum channels for the transmission of 
classical information. Recently, much attention has been devoted to the study of quantum channels with memory (see e.g. |Q]] 
and the references therein). There are evidences that memory effects can enhance the classical capacity in both the cases of 
discrete JUl and continuous variables. Bosonic Gaussian channels provide a fertile benchmark for the study of quantum 
channels with continuous alphabet J4j5|]- Here we consider a model of lossy bosonic memory channel introduced in J6|], based 
on the general scheme proposed in |0](see also (Hi). In this kind of model, memory effects come from the interaction with a 
common environment containing correlations. In spite of the fact that each channel belonging to this family is unitary equivalent 

<— — i to a memoryless one, the presence of energy constraints can break the unitary symmetry, leaving the problem of evaluating 
capacities open (see 01). A specific instance of a channel belonging to that family is obtained by specifying the state of the 
environment. Here we consider a multimode squeezed vacuum with one free parameter expressing the degree of squeezing, this 
parameter in turn determines the amount of memory contained in the channel. 

The aim of this paper is to evaluate the maximum achievable transmission rate using homodyne detection. A similar analysis 
^ ' was already presented in where the maximization was performed over a specific set of encoding/decoding schemes; here 
O" 1 we optimize the rate over a much larger set. That allows to compare explicitly the performance of suitably defined collective 
l— ~ 1 ■ schemes, i.e. using entangled states, over local schemes involving only separable ones. 

| The paper proceeds as follows. In section|ll] we introduce the model. In section [Til] we describe the procedure for encoding 

■ and decoding of classical information. The main results are presented in section HVl where we compare local and collective 
\Q ' encoding/decoding schemes. SectionlVlcontains conclusions and comments. 
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"xf \ II. LOSSY BOSONIC CHANNEL WITH MEMORY 

To define our model we need to introduce, for a given integer n, a set of n input bosonic oscillators, with canonical vari- 
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ables {qkiPk}k=i,...m an d a collection of n ancillary modes, which play the role of the environment, with canonical variables 
{Qk, Pk]k=i,...n- In the following we refer to this set of oscillators as the local modes. All the frequencies are assumed to be 
degenerate and normalized to one. The integer k labels the sequential uses of the channel. At the feth use, the fcth input mode 
is linearly mixed with the fcth environment mode at a beam splitter with given transmissivity 77 (see figure [T). In the Heisenberg 
picture, the channel transforms the input field variables as 



P'k = VvPk + VT^Pk- Kl} 

Memory effects among different uses of the channel are present if the corresponding local modes of the environment are 
correlated. Here we assume the environment to be in a Gaussian state with zero mean, described by the Wigner function 

W= - 1 =exp (-lx T V~ 1 X) 7 (2) 



y/det(V) \ 2 

where X T := (Qi, Q2, ■ ■ ■ Qn, Pi, P2, ■ ■ ■ Pn) indicates the vector of coordinates in the environment phase space. The choice 
of a Gaussian state for the environment makes the channel itself Gaussian. We chose a covariance matrix with the following 
block-diagonal form 
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FIG. 1: A schematic picture of the model of lossy bosonic channel. Each input mode (left-right line), representing one use of the channel, 
interacts with the corresponding environment mode (top-bottom line) through a beam-splitter. To introduce memory effects, environment 
modes are considered in a correlated state. 



This is a bona fide covariance matrix as long as the matrix Q is symmetric and the parameter s is real. It represents a multimode 
squeezed vacuum (see e.g. J^t). The squeezing parameter s (or |s|), determining how strong environment correlations are, can 
be also interpreted as a measure of memory between different uses of the bosonic channel. As to the form of fi, we chose the 
following n x n matrix: 
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The spectrum of this matrix was already presented in e.g. |@], its eigenvalues are 

nj 



A, = 2 cos 



the components of the corresponding eigenvectors being 
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for j = 1, . . .n, 
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for k = 1, 
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According to Williamson theorem (see e.g. 11010 one can always find a set of collective modes, with canonical variables 
{Qj, Pj}j=i.... n , which diagonalize the covariance matrix. In this basis, the Wigner function has the form 
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where X T := (Qi, Pi, ■ ■ ■ Q n , Pn)- An explicit expression for these variables can be obtained from the eigenvectors of the 
matrix f2, yielding to define the following variables: 



Pj ■= Vj, k Pk- 

Analogously, we introduce the following collective variables at the input field: 
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We refer to the set of oscillators with field variables {Qj, Pj}j=i,... n an d {Qj,Pj}j=i,...n as the collective modes. Let us remark 
that the action of the channel in the Heisenberg picture on the collective modes is formally unchanged: 



Qj = Vv Qj + y i - V Qj , 
Pj = VVPj + V 1 -V Pj ■ 



(10) 
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It is worth noticing that the diagonal covariance matrix V is the direct sum of one-mode covariance matrices: 

n 
3=1 
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That implies that the state of the environment is simply separable in the basis of collective modes. It follows from the trans- 
formation (O that the entries of V are the eigenvalues of V. Using the eigenvalues of the matrix il, we obtain the following 
expression for the covariance matrix of the jth collective mode: 
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where Sj := sXj. We observe that each collective oscillator of the environment is in a squeezed state with j dependent squeezing 
parameter. 

To conclude this section, let us come back to the basis of the local modes. For a given integer fco G [1, n], by integrating the 
Wigner function in (0 over the local variables {qk,Pk} for k ^ fco we obtain the Wigner function describing the state of the 
fcoth local mode of the environment. The corresponding local state is thermal-like, with average number of excitations: 



T e g(s, fc ) = i 



3=1 



(13) 



By the symmetries of Vjk and Sj, this can be rewritten as follows: 



T ff(s, fco) 
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j,fco COShs 3 
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Ej=i 1)/2 ^ 2 ,fc coslls J 

Notice that the local temperature T e s(s, fco) is a monotonically increasing function of \s 



if n is even, 
if n is odd. 



(14) 



III. CONTINUOUS VARIABLES ENCODING/DECODING 



Classical information is sent via the bosonic channel by choosing a suitable scheme for encoding (by state preparation of the 
input field) and decoding (by observing the output field) a classical alphabet. 

As to the encoding, we introduce a reference pure state of the input field, described by the Gaussian Wigner function 
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expressed in terms of the local variable x 



. q n ,pi, ■ ■ -p n 



A multivariate Gaussian variable y 



(y q ,i, ■ ■ ■ y q , n ,yp,i, ■ ■ ■ yp,n) T , with probability distribution P(y), zero mean and covariance matrix Y, is encoded in the dis- 
placed state of the input field described by the Wigner function 



y/det((r) 
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Eventually, the state describing the statistical ensemble has the following Wigner function: 
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To avoid infinite energy we introduce a constraint in the maximum number of excitations at the input field per channel use 
(i.e. per mode) in average. Allowing no more than N excitations per mode in average, the constraint can be written in terms of 
the covariance matrices, in natural units: 



-Ltr(a + Y)<N+l 
In 2 



(18) 
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We can summarize the action of the channel on Gaussian states as follows. For a given value of the displacement amplitudes 
(hence for a given letter of the alphabet) the Wigner function at the output field has covariance matrix 

a' y = W + (1 - n)V. (19) 

On the other hand, the output state averaged over the letters of the continuous alphabet has covariance matrix 

a' = r)(a + Y) + (1 - r))V. (20) 

We can equivalently work in the basis of collective modes, in which the covariance matrix of the reference input state is 
denoted as a. Analogously, the zero mean Gaussian variable y :~ (y q .i, ■ ■ ■ Vp.n) 1 , defined by the relations 



Vq.j ~ V j,kVq,k, 

V P ,j ■= v hk y Ptk , 



(21) 



has covariance matrix Y . In the collective basis, the output field is described by the covariance matrices 

a' y =w + (l-r))V (22) 

and 

a' = r)(a + Y) + (1 - V )V. (23) 
It is worth noticing that the form of the energy constraint is preserved, namely 

— tr (a + Y) <N + -. (24) 



2n 

Homodyne detection requires the choice of a compatible set of n quadratures to be measured at the output field. A generic set 
of quadratures is 

{z h | z h =R hk q k +S hk p k } (25) 



for any pair ofnxn matrices R, § satisfying the relations RR T + §§ T = I and RS T - SIR 1 " = O (see e.g. (TOD)- Assuming 
ideal homodyne, the distribution of the stochastic variable z := (z\, . . . z n ) T is Gaussian with zero mean. Its covariance matrix, 
which we denote as Z, can be computed from the output field covariance matrix using the relations in (l25l . 

For given covariance matrices V, a, and for a given set of quadratures to be measured, a classical channel is defined from the 
quantum one. The capacity of the classical channel can be computed as the maximum of the mutual information between the 
output and the input variables: 

I(z;y)=H(z)-H(z\y), (26) 

where H denotes the Shannon entropy. The maximization has to be taken over all possible expressions of the covariance matrix 
Y compatible with the energy constraint. 



IV. ACHIEVABLE BIT RATES WITH HOMODYNE DETECTION 



Let us now consider the form of the Wigner function for the reference input state, determined by the covariance matrix a, and 
the form of the distribution of the classical variable y, determined by the covariance matrix Y . Here we distinguish and compare 
two cases: the first one corresponds to a local encoding in which a and Y are diagonal; the second one is a collective encoding 
in which a and Y are diagonal. From the view point of the local modes, in the local encoding scheme information is always 
carried by simply separable (unentangled) states, while the collective encoding deals with states which are in general entangled. 

Let us first describe the case of local encoding. In the basis of the local modes, the diagonal a can be parameterized as follows 

n 

<7 = 0a fc . (27) 

k=l 

We introduce the real parameters {r k }, such that 
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Analogously, the covariance matrix Y reads 



(29) 



k=l 



and is parameterized by the positive parameters {c q , k ,c p . k }: 



1 ( c q , k 



c Pifc 

As to the decoding part, the natural choice is to measure the local quadratures 

z k := cos k q k + sm6 k p k . 



(30) 



(31) 



The covariance matrix of the output variable z := [z\, . . . z n ) T can be easily computed. It is diagonal as well, with the 
variances: 



cos k 



V [ ^tt + c q ,k) + (1 - V) [T cS (s,k) + 1 



sin 2 k 



+ c P ,k + (1 - r?) T eff (s, k) + - 
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The output variable has conditioned covariance matrix, denoted Z y , which is diagonal with entries 



(4)v = cos 2 k 



sin k 



+ (!-»?) [T eS ( S> k) + - 



(32) 
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Since all the probability distributions are Gaussian, the mutual information can be easily written in terms of the covariance 
matrices: 



1 1 1 ™ 

I(z; y) = - log 2 [det(Z)] - - log 2 [det(Z„)] = - log 2 



fc=i 



(4)v 



(34) 



The maximization of the mutual information is over the 3n parameters r k , c qtk , c Ptk satisfying the constraints 

1 " 1 

^- (coshrfc + c q>k + c Pifc ) = N k + -, 
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(35) 
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and over the n angles k , which determine the chosen quadratures. Concerning the choice of the optimal quadratures, it is 
immediate to recognize that the maximum is reached for sin 6 k = if r k > 0, and for cos k = if r k < 0, while the value of 
9 k is irrelevant if r k = 0. Maximizing over the remaining variables and applying the additivity of mutual information, we obtain 
the following expression for the maximum mutual information with local encoding/decoding scheme: 



-Fi oc = max 
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2N k + 1 - cosh r° pt 
e~ r T +u(2T eS {s,k) + l) 
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where 



e fc 



i/(2T cff (s,fc) + l) 



VI + v{2T cS (s, k) + 1) (^(2T cff (s, k) + 1) + 4Ar fe + 2) - 1 



(37) 



(38) 



and we have defined v := (1 — r))/r]. 

Let us now consider the case of collective encoding. In the basis of the collective modes, the diagonal a can be parameterized 
as follows 



* = °3 

3 = 1 



(39) 
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with 



e r = 
e~ f 



Analogously, the covariance matrix Y is 



with 



Y, 
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As to the decoding part, the natural choice is to measure the collective quadratures defined as 

Zj = cos 6 jqj + sin 6 jpj . 



(40) 



(41) 



(42) 
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The output variable z := (z%, . . . z n ) T has diagonal covariance matrix Z and conditional covariance matrix Z y . The correspond- 
ing variances read 



(z 2 ) = cos 2 I 



and 



<«i)» =cos 2 l 



sin 



sin 2 0o 



I ^-+cp,j ] +(l-»7)- 



^7^- + (1 - ? 7)^- 



(44) 
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It is easy to recognize that the maximum of the mutual information is reached for sin 8j — if Sj > 0, and for cos 8j = 
if Sj < 0, finally if Sj = we can argue as in the case of local encoding. By maximizing over the remaining in variables {fj, 
c q j, c p j} under the constraints 



1 - 1 
— 2^ (cosh ^ + cgj + £pj) = Nj + - 

i - 



3=1 



(46) 



(47) 



3=1 



and applying the additivity of the mutual information, we obtain the following expression for the maximum rate of transmission 
with encoding/decoding in collective variables: 



F a 



max < - log 2 



1 + 2- 



2Ni 



1 — cosh f ° pt 
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< JV 



3=1 



The optimal value f° pt is determined by the equation 



1 + ve~\ s i\ (ve~\ s i \ + AN 3 + 2) - 1 



(48) 



(49) 



The maximum mutual informations per channel use, respectively defined as F\ oc /n and F co \\/n, are plotted in figure [2] as 
function of number of uses of the channel, for several values of the memory parameter. The case n = 1 or s = correspond to 
the memoryless channel. Notice the different behavior of the rates as the number of uses (or the value of the memory parameter) 
increases. The same quantities are plotted together in figure 13 for n = 10 uses of the channel, as function of the memory 
parameter, for several values of the beam splitter transmissivity. Notice the different behavior of the rates corresponding to local 
(dotted line) and collective (solid line) encoding/decoding scheme. 

The expression of the maximal mutual information allows to make comparison with the memoryless case (see also the plot 
in figure [3]). First of all, let us recall that the memoryless channel is recovered for s = 0; our results show that the maximum 
mutual information is a monotonically increasing function of \s\ in the case of collective encoding/decoding scheme, while it is 
monotonically decreasing in the case of the local one. This is enough to conclude that, in the presence of memory, collective 
encoding/decoding is the optimal choice. Secondly, let us notice that the two schemes are equivalent in the memoryless case, i.e. 
s = orn = 1; that leads to conclude that the presence of memory enhances the bit rate as long as collective encoding/decoding 
is allowed. 
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FIG. 2: The histograms show the maximum achievable rates per channel use with collective (a) and local (b) encoding/decoding scheme, as 
function of the number of channel uses. For each value of n, the histograms correspond to different values of the memory parameter, from the 
left to the right s = 0, s = 1, s = 2, s = 3. The maximum average excitations per mode is N = 8, the transmissivity r\ — 0.7. 




FIG. 3: The plot shows the achievable rates per channel use with collective (solid line) and local (dotted line) encoding/decoding scheme, as 
function of the memory parameter, for n = 10 uses of the channel. The curves correspond to different values of the transmissivity, from the 
bottom to the top r\ varies from r\ — 0.1 to r\ = 0.9 by steps of 0.1. The maximum average excitations per mode is N — 8, which corresponds 
to the asymptotic value (plotted in dashed line) of the rate with collective scheme (see equation d51t): lim,,^^ F co \\/n = log 2 (2iV + 1) ~ 
4.0875. 



V. CONCLUSION 

We have discussed the model of lossy bosonic Gaussian channel with memory presented in @], an instance of a general 
class of bosonic memory channel introduced in 10]. Our model is characterized by a single parameter which determines the 
amount of memory contained in the channel. By considering the case of homodyne detection, we have compared two different 
encoding/decoding schemes: a local one in which classical information is encoded and carried by unentangled (simply separable) 
states; and a collective scheme based on states which are in general entangled among different uses of the channel. We have 
computed the maximum achievable bit rate using both the schemes, assuming energy constrains at the input field, for an arbitrary 
number of channel uses. Our results lead to conclude that, as the memory parameter increases, the collective scheme becomes 
more and more efficient than the local one. It is hence interesting to analyze the asymptotic limit of the rates for |s| — ► oo, 
corresponding to infinite memory. From the expressions in d37| i and d48l l we obtain 

Um — = (50) 

and 

Um =log 2 (27Y + l). (51) 
s— too n 
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Notice that the latter expression holds independently of n and for any value of 77 ^ 0,1. Interestingly enough, the asymptotic 
value in ( |5TT i equals the maximum achievable rate of the (memoryless) noiseless channel (r/ = 1) with homodyne detection (see 
J3l). In this sense, by using local encoding/decoding one sees infinite noise in the channel with high degree of memory, while 
by using a collective encoding/decoding one can completely avoid it. 
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